OTIC  FILE  COPY 


*  * 


RADC-TR-90-1 55 
Final  Technical  Report 
July  1990 


RESONANCE  REGION  SCATTERING 


Syracuse  University 
Xiaopu  Yang  and  Ercument  Arvas 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


Rome  Air  Development  Center 
Air  Force  Systems  Command 
Griff iss  Air  Force  Base,  NY  13441-5700 

vi)  Ja  ^ 


RADC-TR-90-155  has  been  reviewed  and  is  approved  for  publication. 


EDWARD  M.  STARCZEWSKI 


Project  Engineer 


APPROVED: 


FRED  3.  DEMMA,  Actg  Director 
Directorate  of  Surveillance 


FOR  THE  COMMANDER: 


Directorate  of  Plans  &  Programs 


DESTRUCTION  NOTICE  -  For  classified  documents,  follow  the  procedures  in  DOD 
5200.22-M,  Industrial  Security  Manual,  or  DOD  5200.1-R,  Information  Security  Program 
Regulation.  For  unclassified,  limited  documents,  destroy  by  any  method  that  will 
prevent  disclosure  of  contents  or  reconstruction  of  the  document. 

If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization,  please  notify  RADC 
(OCTrvU,  Griffiss  AFB  NY  13441-5700.  This  will  assist  us  in  maintaining  a  current 
mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  requires  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 

Form  Appro  rod 

OPMNo.  0704-019$ 

’  AGENCY  USE  ONLY  (iMn  a  mot  !  2.  REPORT  DATE  3.  REPORT TYPEANO  OATES  COVERED 

3uly  1990  FinaJ  jun  g9  _  Feb  90 

4.  TITLE  AND  SUBTITLE 

RESONANCE  REGION  SCATTERING 

S.  FUNDING  NUMBERS 

C  -  F30602-88-D-0027 

PE  -  6374 ID 

PR  -  3640 

TA  -  04 

WU.-  PI 

VAUTHOR(S) 

Xiaopu  Yang  and  Ercument  Arvas 

7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  ADDRESSEES) 

Syracuse  University 

Dept  of  Electrical  and  Computer  Engineering 

Syracuse  NY  13244 

1.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

N/A 

9.  SPONSOWNG/MONTTORINO  AGENCY  NAMEESi  AND  ADDRESSEES) 

Rome  Air  Development  Center  (OCTM) 

Griffiss  AFB  NY  13441-5700 

10  SPONSORINGMONITORING  AGENCY 
REPORT  NUMBER 

RADC-TR-90-155 

11 .  supplementary  notes 

RADC  Project  Engineer:  Edward  M.  Starczewski/OCTM/(315)  330-4431 

12*.  DISTRIUUTION/AVAILABICrTY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

12B.  DISTRIBUTION  CODE 

1 1  ABSTRACT  (tMxmjm  300  worm} 

Using  frequency  derivative  information,  techniques  were  developed  to  predict  the  scattering 
from  conductor  and  dielectric  2-D  targets  over  a  frequency  band.  Instead  of  computing  a 
method  of  moments  solution  using  a  point  by  point  approach,  a  rational  function  model  is  used 
to  approximate  the  current  as  a  function  of  frequency.  The  model  coefficients  are  computed 
using  frequency  derivative  information  at  one  frequency  within  the  band  to  improve  efficiency 
of  conventional  moment  method  approaches.  Results  show  that  the  model  based  approach  gives 
excellent  results  over  a  limited  frequency  band,  and  is  much  more  efficient  than  the 
conventional  point  by  point  approach.  //,.  ,  S) 

14.  SUBJECT  TERMS 

'f  Spectral  Estimation  ,  Method  of  Moments,  i 

Frequency  Derivatives  ,  Radar  Cross  Section*  J 

IV  NUMBER  OF  PAMS 

36 

IV  PRICE  COOt 

— — WI^W 

<0  LIMITATION  OF  ABSTRACT 

UL 

NSN  7540*01 -2W.5500 


NW  »»  *W  V*  *»>% 

hm\ 


Use  of  Frequency-Derivative  Information 
to  Compute  EM-Scattering  from  Multiple 
Perfectly  Conducting  and  Dielectric 
Cylinders  of  Arbitrary  Cross-Section 
over  a  Frequency  Band 

r  y 

by  Xiaopu  Yang  :  ‘r 

Ercument  Arvas 


I .  Introduction 

The  use  of  frequency-derivative  information  incorporated 
with  model-based  parameter  estimation  (MBPE)  to  obtain 
electromagnetic  response  of  an  antenna,  propagation  path,  or  a 
scatterer  over  a  spectrum  of  frequencies  has  been  studied  by  a 
few  researchers  [1] .  In  this  report,  we  present  some 
interesting  results  obtained  when  the  theory  is  applied  to  E- 
field  solution  of  TM-Scattering  from  multiple  perfectly 
conducting  and  lossy  dielectric  cylinders  of  arbitrary  cross¬ 


section  . 


The  moment  method  solution  of  above  mentioned  problem  has 
been  formulated  in  [2] .  When  the  solution  over  a  frequency  band 
are  required,  the  method  used  most  often  is  the  one  that  finds 
the  solution  at  a  set  of  discrete  frequencies,  and  then 
interpolates  them  over  the  frequency  band.  The  cost  of 
computation  for  this  conventional  method  can  be  very  high  when 
the  number  of  discrete  frequencies  is  high  as  it  is  usually  the 
case  to  obtain  esthetic  appearance  of  low-order  interpolatory 
curve.  In  contrast,  MBPE  approach  exploits  the  underlying 
physics  of  the  phenomenon  being  modeled  so  that  unnecessary 
computation  is  avoided.  Using  frequency-derivative  information 
further  reduces  the  cost  of  computation  since  as  far  as 
determination  of  model  coefficients  is  concerned,  frequency- 
derivative  samples  bear  the  same  amount  of  information  as 
frequency  samples,  but  the  former  are  much  easier  to  obtain. 
Moreover,  as  demonstrated  by  numerical  results  presented  in 
this  report,  using  frequency-derivative  information  has  other 
advantages,  such  as  better  local  approximation  and  easier 
implementation . 


II.  Model-Based  Parameter  Estimation 


A.  A  Rational  Function  Representation 
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The  MBPE  approach  utilizes  the  mathematical  property  that 
the  frequency  response  is  well  approximated  by  a  rational 
function  of  form 

F  (s)  *  N(s)/D(s)  (1) 

P  Q 

-  [I  NiS1]/!!  DiS1], 

i=0  i=0 

which  in  turn  is  originated  from  the  Singularity-Expansion 
Method  [3]  in  using  a  pole  series  in  complex  frequency  s*0+jw 
to  model  the  frequency  response,  i.e. 


M  N 

F(s)=  Fp  (s)  +Fnp  (s)  =  I  R«/(s-s0)  +  I  CpsP  (2) 

a-1  P=-L 

where  Fnp(s)  accounts  for  the  nonpole  portion  of  F(s),  and  in 
Fp(s)  the  parameters  Rq  and  sa  are  the  residues  and 
corresponding  poles.  While  the  number  of  poles  is  theoretically 
infinite  to  cover  the  entire  frequency  band,  quite  often  an 
approximation  to  F(s)  over  a  limited  frequency  band  is  of 
interest  so  that  the  number  of  poles  used  in  the  model  is 
small,  and  problem  of  ill  conditioning  can  be  avoided.  More 
models  can  be  used  to  cover  a  wider  frequency  band.  In  our 
results,  surface  currents  defined  in  the  method  of  moments  are 
modeled  by  a  rational  function  of  frequency.  For  simplicity, 
unless  we  mention  otherwise,  the  current  is  modeled  by 
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1(a) 


Ra/ (S-Sa) +Cq, 


(3) 


P 

V 

4. 

a-1 

which  is  equivalent  to 

I  (s)  -  N(s)/D(s) 

P  P 

-  NiS1]/!!  Dts1],  (4) 

i=0  i=0 

in  which  there  are  2P+1  coefficients  N^'s  and  Dj/s  to  be 
determined  (Do=l) • 

B.  Computing  Model  Coefficients  Using  Frequency-Derivative 
Samples 

The  coefficients  N-l's  and  Dj/s  of  Eq.  (4)  can  be  found  when 
frequency-derivative  samples  at  some  frequency  (s^so^jwo)  are 
available,  where  w0  is  in  the  frequency  band  of  interest. 
Starting  with  Eq. (4)  and  differentiating  t  times  with  respect 
to  s,  there  results  the  following 

I  ( s ) D  ( s ) =N ( s ) 

I '  (s)D(s) +1  (s)D'  ( s )  =N '  (s) 

I”  (s)  D  (s)  +21'  (s)D'  (s)  +1  (s)D"  (s)=N"  (s)  (5) 
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I  (t)  (s) D  <s)  +tl  <t_1)D'  (s)  +  ...  +Ct,  t_mI  (m)  (s)  D  (t-m)  + 

+. . .+1 <s)D(t)  (s)=N<fc)  (s) 

where  Cr, a  =  r ! / [s !  (r-s)  ! ]  is  the  binomial  coefficient.  The 
system  of  t+1  equations  in  Eq. (5)  provides  the  information  from 
which  the  model  coefficients  can  be  found  if  t+l£2P+l.  With  no 
loss  of  generality,  Eq. (5)  can  be  simplified  by  setting  so=0 
with  s  representing  the  frequency  variation  about  s0-  We  then 
obtain  the  following  matrix  representation  for  the  unknown 
coefficients 


-1  0...0  0  0...  0- 

r  N0  - 

r  *o  1 

0  1  ...  0  -I0  o  ...  0 

Ml 

II 

0  0  ...  0  -Ii  -Iq  ...  o 

n2 

I2 

0  0  ...  1  -Ip - 1  - Ip  — 2  *  *  .  -Io 

Np 

Ip 

0  0  ...  0  -Ip  -Ip-i ... -Ii 

Di 

= 

Ip+l 

0  ...  0  -Ip+i  . . .  -I2 

d2 

Ip  +  2 

0  ...  0  -I2P-1  • • *  -ip 

Dp 

I2p 

where  Im  =  (1/m! ) I <m) (s0) . 

C.  Getting  Frequency-Derivative  Information  from  more  than  one 
Points . 

For  purpose  of  comparison,  numerical  results  are  also 
presented  of  method  using  frequency-derivative  information  from 
more  than  one  frequency  points.  Specifically,  in  two  points 
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frequency  derivative  method,  frequency  samples  and  frequency 
derivative  samples  at  s0±s  are  used  to  find  model  coefficients. 
And  in  three  points  frequency  derivative  method,  the 
informations  are  taken  at  So  and  Sq±s. 


By  following  the  concept  in  Sec  II. B,  we  can  obtain 
similar  matrix  equations  for  the  model  coefficients.  For 
example,  when  frequency  samples  up  to  the  second  order 
derivatives  at  two  points  are  used,  the  matrix  equation  which 
determines  coefficients  in  model 


I(s)  =  (No+NxS+NzS2) / (1+D1s+D2a2+D3s3)  (7) 


is 


1 

-s 

s2 

I0 (-3) S 

-I0  (-3)  S2 

0 

1 

-2s 

II  ( ~  S )  3- 

"II  (-3)  S2  + 

-I0  (-3) 

+2I0  (-s) S 

0 

0 

1 

I2  (-S)  3- 

-I2  (-3)  S2  + 

-II  (-3) 

+  21],  (  —  S  )  s- 
-I0  (-3) 

1 

s 

s2 

-lo  (3)  3 

-I0 (3) S2 

0 

1 

2s 

-Ii<s)s- 

-Ii  (3)  S2- 

-I0  (3) 

-2I0  (3) s 

n 

0 

1 

-I2  (s) s- 

-I2  (3)  s2- 

-II  (3) 

-2IX  (s) s- 
-Io  (3) 

Io(-s)s3 

■  No  - 

-Io(-s)' 

l!(-S)S3- 
-3I0 ( - s ) s2 

Nx 

lx (-3) 

I2(-s)s3- 
-3Ii (-s) s2+ 
+  3I0 (-s) s 

n2 

I2  (-s) 

-I0 (s) s3 

Di 

Io  (3) 

-Ii  (s) s3- 
-3I0  (s) s2 

d2 

Il(3) 

-I2 (s) s3- 
-3Ii  (s) s2- 

d3 

1 2  (3) 

-3I0 (s) s 


where  Im 


(l/m!)I(m)  ( s o )  •  Similarly,  we  can  find  the  matrix 


equation  for  coefficients  of  model 
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I(s)  -  (No+N13+N2a2+N3s3+N4s4) / (1+D1s+D2s2+D3s3+D4s4)  (9) 

using  information  up  to  the  second  order  derivatives  from  three 
points . 


III.  Computing  Frequency  Derivatives  in  A  Moment-Method 

Model 

A.  The  Basic  Idea 

Overall  computation  in  obtaining  frequency-derivative 
information  needed  in  Eq. (6)  can  be  dramat ically  reduced  from  a 
moment -method  model. 

On  writing  the  moment -method  equations,  we  have 


N 

I  Zijlj  =  V±  (10) 

j  =  l 


where  Z,  I,  and  V  are  the  impedance  matrix  and  the  current  and 
voltage  vectors.  A  solution  for  the  current  is 


N 

I±  =  I  YijVj  (11) 

j=l 


where  matrix  Y  is  the  inverse  of  Z.  If  we  differentiate  Eq.  (10) 
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with  respect  to  frequency,  we  obtain 


N 

I  (ZijI' j+Z' ijlj)  =  V' j  (12) 

j-1 

from  which  results  the  differentiated  current 

N  N 

I' ±  =  I  YijfV' j-I  Z' j)cIk] .  (13) 

j=l  k=l 

It  is  seen  that  the  differentiated  current  requires  an 
additional  number  of  computations  beyond  those  needed  for  the 
current  proportional  to  N2  rather  than  N3  that  would  be 
required  to  obtain  another  frequency  sample. 

In  general,  the  n'th  frequency  derivative  of  the  current 
is  given  by 


N  N  N 

Ii(n)  =  I  Yi-jCVj  <">-!  Cn,m(Z  Zjk(m)I)c(n-m>)  ]  .  (14) 

j=l  m=l  k=l 

B.  Implementation  in  the  E-field  solution  of  TM-scattering 

Following  the  development  in  [2],  we  have  the  moment 
equation 

Z  I  =  V  (15) 

where 
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z 


(16) 


Zee  ZCD  ZCM  1 
Zdc  Zdd  ZDm  j 
0  Zmd  Z mm  , 

I  =  [  Ic  ID  K°  ]T, 


(17) 


and 


V  =  (  V°  CD  0  ]  T  (18) 

are  the  equivalent  impedance  matrix  and  current  and  voltage 
vectors,  and  they  are  all  functions  of  frequency  either  through 
k0  (=wV|i£  q)  or  k  (=vr/px~)  . 

It  is  seen  from  Eq. (13)  that  in  addition  to  the  impedance 
matrix,  the  differentiated  impedance  matrices  are  needed  in 
obtaining  the  differentiated  currents.  Basically,  there  are  two 
integrations  that  are  numerically  evaluated  in  forming  the 
impedance  matrix.  They  are 


Ai 


r 

Ho2  (klf-f '  I  )dl'  , 


(19) 


and 


A2  =  -jk/4 


jHl- 


<kl f-f'  I  ) 


nid. (f-f' ) / i dl' . 


(20) 


Similarly,  the  differentiated  impedance  matrices  involve 
derivatives  of  those  integrals.  By  taking  derivatives  of 
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gq  (19)  and  (20)  with  respect  to  k  and  using  properties  of 
Hankel  function 

H0<2>'(x)  -  -H1<2)(x)/  (21) 

and 

Hx(2),(x)  =  H0 <2)  (x) -Hi  (2)  (x) /x,  (22) 

we  obtain 


dLAi/dk  =  A1/k+(k^/4)|if-f'  lHl(2>  (klf-f'  l)dl',  (23) 

AX"  =  -A!/k2+Ai2/k+(k^/4)  lj-f'  I2Ho<2)  Odl',  (24) 

« 


Ax<3>  -  3A1/k3-3A1' /k2+2A1”/k- 

-(k*J/4)  \f-f'  I  3HX  <2)  ()dl'  ,  (25) 

Ai(4)  =  -9Ai/k4+9Ai' /k3-5A1"/k2+2A1 (3) /k- 

-(k^/4)|l f-f  i4H0(2)  <)dl',  (26) 


Ai<5) 


45A1/ks-45A1' /k4+24A1"/k3-9A1 <3) /k2+3A! <4) /k+ 
+  (k^/4)  jlf-/'  f  5Hi  (2)  Odl' ,  (27) 
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A-,<6>  -  -225A1/k6+225A1' /k5-127A1"/k4+42A1 <3) /k3- 

-12AX  <4)  /k2  +  3Ai  (5)  /k+ 

+  <k»f/4)  jlf-f'  i6H0(2>  Odl',  C 


1575A1/k7-1575A1/ /k6+810Ai"/k5-285A1 <3> /k4  + 
+78AX <4) /k3-18Ai (5) /k2+4Ax  <6) /k- 


-  (klj/4)  lf-f'  l7H1  <2)  0  dl' 


-HO25A1/k0  +  llO25Ai'  /k7-5625Aj.  ”/k6+ 
2950Ai (3) /k5-519Aj/4' /k4+114Ai (5) /k3- 
-22AX <6) /k2+4Ax  <7'/k- 


-  (k1£/ 4)  !f-f'  I  8Hq(2)  Odl'  , 


dA2/dk  »  -jk/4  Ho<2) (klf-f' l)n±d. (f-f' )dl' 


A2'  +  (jk/4)  Kx  <2)  ()  [r»id.  (f-f'  )  ]  If-/'  Idl' 

J 


A2<3)  =  -A2'/k2+A2"/k+ 

+ ( jk/4)  Ho<2)  0  [nid.(f-f')]  lf-f'  l2dl',  (33) 


A2<4)  ■  3A2' /k3-3A2"/k2+2A2 <3) /k- 

-(jk/4)  HX<2>  ()  [ru'Mf-f')]  lf-f'  l3dl',  (34) 

A2<5>  .  -9A2' /k4+9A2"/k3-5A2 <3) /k2+2A2 (4> /k- 

-  (jk/4)  [h0  (2)  0  [nid.  (f-f'  )  1  If-/'  I  4dl'  ,  (35) 
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A2<6)  -  45A2' /k5-45A2"/k4+24A2 <3> /k3-9A2 <4) /k2+ 

+3A2 <5) /k+ 

+  <jk/4)  HX  <2)  0  [nid.  (f-f  )  ]  \?-f  l  5dl'  ,  (36) 


A2<7)  *  -225A2' /k6+225A2'7k5-117A2 <3> /k4+42A2 <4) /k3- 
-12A2<5>/k2+3A2<6>/k- 

”  ( jk/ 4)  j^Ho  <2)  ()  [n±d.  {f-f  )  ]  I  f-f'  l6dl',  (37) 

A2<8)  =  1575A,'  /k7-1575A2"/k6  +  910A2 <3> /k5-285A2 <4> /k4  + 
+78A2 <5) /k3-18A2 '6) /k2+4A2 <7> /k+ 

+  ( jk/4) |h! <2)  ()  [nid,  if-f')]  I 7dl'  .  (38) 

Using  Eq.  (23) -  (38) ,  the  differentiated  impedance  matrices  take 
the  form 


z'cc 

2'  CD 

Z'  CM 

z '  ■ 

Z'dc 

2  '  DD 

Z  '  DM 

L  0 

2'  MD€ 

ZW 

r 

r  7 " 

^  CC 

2"  CD 

Z  ”  CM 

- 

z"  = 

7  M 

L  DC 

Z"dd 

.  Z  "dm 

L  0 

Z’W 

2  Z"mm£ 

(39) 


(40) 


etc.,  where  E  is  dielectric  constant  of  the  dielectric  body 


IV.  Sample  Numerical  Results 

Included  in  this  section  are  sample  results  obtained  from 
a  Fortran  computer  code  which  implements  the  theory  presented 
above.  Specifically,  the  frequency  derivative  method  (FD)  is 
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compared  with  conventional  point  by  point  method  (CV)  in 
obtaining  the  backscattering  cross-section  (BSCS)  of  various 
structures  over  a  frequency  band. 

Fig. (1)  shows  the  BSCS  of  a  conducting  cylinder  over  a 
frequency  band  from  k=0.5ko  to  k=2k0,  where  k  is  wave  constant, 
and  k0  corresponds  to  the  frequency  at  which  frequency 
derivative  samples  are  taken.  There  are  two  curves  for  FD 
method.  One  uses  frequency  derivative  samples  up  to  4th  order, 
another  up  to  6th  order.  It  is  seen  that  the  curve  for  6th 
order  FD  method  is  in  excellent  agreement  with  that  of  CV  over 
entire  frequency  band  of  interest,  while  the  curve  for  4th 
order  FD  method  starts  deviating  around  k=1.7k0. 

Fig. (2)  shows  the  BSCS  of  a  dielectric  cylinder  (£r*4)  of 
rectangular  cross-section.  The  dimensions  of  the  cylinder  and 
the  number  of  sections  approximating  the  surface  in  applying 
moment  method  are  given  in  the  inset  of  the  figure.  The 
observation  is  similar  to  that  of  Fig. (1) . 

Consider  a  combination  of  conductor  and  dielectric 
cylinders  as  shown  in  Fig. (3) .  The  curves  of  BSCS  for  FD  method 
are  in  fairly  good  agreement  with  that  of  CV  method  within  the 
neighborhood  of  k0  (k=0.8ko  to  k=1.2ko) •  The  6th  order  FD 
method  gives  better  approximation  than  the  4th  order  FD  within 
frequency  band  from  k=0.68ko  to  k=1.28k0. 
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An  interesting  example  is  the  BSCS  for  an  empty  circular 
cylinder  with  an  opening  00*5°  .  The  curve  for  CV  method  is 
presented  in  [4] .  There  are  two  resonances  within  frequency 
band  from  ka=4  to  ka=6,  where  a  is  the  radius  of  the  cylinder. 
Fig. (4)  shows  the  curves  for  FD  method  with  frequency 
derivative  samples  taken  at  ka=5  (k0a«5)  along  with  that  for  CV 
method.  The  FD  method  of  4th,  6th,  and  8th  order  are  used.  All 
three  curves  of  FD  method  locate  the  resonance  around  ka=5.16. 
Notice  that  the  higher  the  order  of  FD  method,  the  better  its 
curve  agrees  with  CV  curve.  In  particular,  the  4th  order  curve 
misses  the  resonance  at  ka=5.55  completely,  while  the  6th  order 
curve  gives  some  indication  of  resonance  there,  and  the  8th 
order  curve  locates  the  resonance.  Next,  we  take  frequency 
derivative  samples  at  ka»5 . 3  (k0a*5.3).  Then  as  shown  in 
Fig. (5) ,  all  three  curve  of  FD  are  able  to  pick  up  the  existing 
two  resonances.  However,  only  8th  order  FD  curve  is  in 
excellent  agreement  with  the  CV  curve.  Finally,  several  models 
are  used  to  cover  a  frequency  band  from  ka=0.6  to  ka=6.6. 
Specifically,  entire  frequency  section  is  divided  into 
subsections  of  equal  size,  and  the  solution  in  each  section  is 
modeled  by  a  rational  function  whose  coefficients  are 
determined  using  6th  order  FD  method.  For  instance,  frequency 
derivative  samples  taken  at  ka=0 . 8  are  used  to  cover  the  region 
from  ka®0 . 6  to  ka=1.0,  samples  taken  at  ka=1.2  are  used  to 
cover  ka*l . 0  to  ka=1.4,  and  so  on.  The  combined  result  of  those 
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models  are  shown  in  Fig. (6)  which  gives  the  same  locations  of 
resonance  as  given  in  Fig. 8(b)  of  [4]. 

Fig. (7)  shows  the  result  using  two  points  FD  method  for 
the  same  structure  used  by  Fig.  (4) -(6).  Informations  used  to 
determine  model  coefficients  are  taken  at  ka  ™  k0a±5ka,  where 
k0a  *  5.0.  Curves  corresponding  to  different  Ska  are  shown.  It 
is  shown  that  the  curve  varies  with  Ska.  As  far  as  frequency 
band  from  ka=4 .  to  ka=6 .  is  concerned,  except  the  curve  with 
Ska**0.8,  all  three  other  curves  indicate  both  resonances  at 
ka~5.16  and  ka=5.55. 

Fig. (8)  shows  the  result  using  three  points  FD  method.  The 
informations  are  taken  at  ka  =  koa,  and  ka  =  k0a±5ka,  where  kQa 
is  again  5.0.  Among  curves  corresponding  to  four  different  Ska, 
the  one  with  Ska  »  0.3  gives  the  best  result.  In  addition  to 
existing  resonances,  other  curves  show  signs  of  resonances  that 
are  not  existing. 

Finally,  Fig. (9)  and  Fig. (10)  demonstrate  the  effect  of 
using  models  other  than  those  given  by  Eq.  (3)  and  (4) .  What 
shown  are  results  of  6th  order  FD  method  applied  to  previous 
examples  with  different  models.  Models  are  different  in  orders 
of  polynomials  in  its  numerator  and  denominator.  In  both 
figures,  NN  is  the  order  of  numerator,  and  ND  is  the  order  of 
denomenator.  It  is  seen  that  curves  corresponding  to  different 
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models  are  different,  but  they  all  agree  with  each  other 
locally  around  the  point  where  informations  are  taken. 


V.  Summary  and  Conclusions 

The  problem  considered  here  is  that  of  computing  the 
electromagnetic  parameters  over  a  frequency  band.  Instead  of 
using  a  point  by  point  approach,  a  rational  function  model  is 
used  to  approximate  the  parameter  of  interest  as  a  function  of 
frequency.  The  model  coefficients  are  computed  using  frequency 
derivative  information  at  one  frequency  within  the  band  to 
improve  efficiency.  Sample  results  are  presented  which 
illustrate  that  the  model  based  approach  gives  excellent  result 
over  a  limited  frequency  band,  and  is  much  more  efficient  than 
conventional  point  by  point  approach. 

We  have  the  following  conclusions : 

First,  around  the  point  where  frequency  sample  and  its 
derivatives  are  taken,  there  is  always  a  frequency  band  in 
which  the  FD  method  gives  good  approximation  to  the  true  value 
of  the  frequency  response.  The  higher  the  order  of  derivatives 
used,  the  larger  the  size  of  this  frequency  band.  When  the  same 
order  of  derivatives  is  used,  the  size  of  this  band  still 
varies  from  problem  to  problem.  It  depends  on  complexity  of  the 
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function  to  be  modeled,  i.e.  the  difference  between  the  order 
of  the  assumed  rational  function  model  and  the  order  of 
rational  function  with  complete  fitting  to  the  function. 

Secondly,  one  point  FD  is  preferred  to  multi-point  FD  for 
the  following  reasons .  One  point  FD  is  simpler  to  implement  in 
that  it  involves  a  matrix  equation,  Eq.(6),  which  can  be 
constructed  easily,  and  there  is  no  need  to  select  separation 
of  points  as  required  in  multi-points  FD;  One  point  FD  gives 
better  approximation  locally,  and  separate  models  can  be  used 
when  a  broad  frequency  band  is  of  interest. 

Finally,  when  the  sum  of  orders  in  numerator  and 
denomenator  polynomials  is  fixed,  different  combination  of 
orders  generally  yields  different  result.  The  selection  of  the 
best  model  requires  prior  knowledge  of  the  function  to  be 
modeled.  When  such  knowledge  is  not  available,  model  of  Eq.  (4) 
is  recommended  for  its  simplicity. 
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